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Numerical Results for Transport Coefficients of Quark Gluon Plasma 
with Iwasaki's Improved Action 
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b Faculty of Education, Yamamagata University, Yamagata 

Numerical results for the transport coefficients of quark gluon plasma are calculated by lattice simulation of 
SU(3) pure gauge model. The bulk viscosity is consistent with zero. The shear viscosity is finite and increases 
with temperature T roughly as T 3 , and around the finite temperature transition points, it is slightly smaller than 
the typical hadron masses. 



1. Introduction 

In the phenomenological study of quark gluon 
plasma(QGP), when its bulk properties are con- 
cerned, the system of quarks and gluons is usually 
treated as gas or liquid. Then the fundamental 
parameters of QGP such as transport coefficients, 
are very important information. The aim of this 
work is to calculate them from the fundamental 
theory of QCD. 

The calculation of transport coefficients is for- 
mulated in the framework of linear response the- 
ory of Kubo|l]-[5|]. They are expressed by the 
space time integral of retarded Green's function of 
energy momentum tensors at finite temperature. 
The shear viscosity rj is expressed as follows. 
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Similarly the bulk viscosity £ and heat conduc- 
tivity x are expressed in terms of the retarded 
Green's function of Tn and T41 components of en- 
ergy momentum tensor. In the pure gauge model 
the energy momentum tensors are written by the 
the field strength tensor. 

The direct calculation of the retarded 
Green's function at finite temperature is very dif- 
ficult. Then the shortcut is to calculate Matsub- 



ara Green's function(Gyj) and then by the ana- 
lytic continuation, we obtain the retarded Green's 
function at finite temperature. The analytic con- 
tinuation is carried out by the use of the fact 
that the spectral function of Fourier transform of 
the both Green's functions is the same. For the 
spectral function we use the following simplest 
ansatzB. 
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where 7 partially represents the effects of the in- 
teractions and is related to the imaginary part 
of self energy. Under this ansatz, the transport 
coefficients are calculated as, 



a x a s = 2A 
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where a means 77, |?7 + ( and \ ' T. 

We notice that if m = or 7 = 0, trans- 
port coefficient becomes zero. In order to deter- 
mine these parameters, at least three independent 
data points in Gp are required in the tempera- 
ture direction. In the pioneering work of Karsch 
and Wyld||, they performed lattice QCD calcu- 
lations on 8 3 x 4 lattice and the resolution was 
not enough to determine these three parameters 
independently. 
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2. Numerical Results 

We carry out our simulation on 16 3 x 8 lat- 
tice. We have started the the simulation from 
U(l) gauge theory. It is found that the fluc- 
tuation of Gp is very large that it need about 
10 6 data for the determination of Gp in the de- 
confined phase. Further in the confined phase, 
the fluctuations become still larger and we could 
not obtain the Green's functions even with the 
~ 1.5 x 10 6 data|6|. Similar results are obtained 
in the case of SU(2) gauge theory. Then in the 
calculation of SU(3) gauge theory, we carry out 
our simulation only in the deconfined phase. 

It is also found that the fluctuation of Gp be- 
comes larger as we proceed to U(l) — > SU{2) — > 
SU(3). Then it is a very important problem to 
reduce the fluctuation of Gp. We find that by us- 
ing the Iwasaki's improved action, the fluctuation 
is much reduced as shown in Fig.l. Then in the 



timated by the Jackknife method. The results 
for them are shown in Fig. 2. It is found that the 
bulk viscosity is zero within errors. This result 



Figure 2. Shear and bulk viscosity xa 3 



Figure 1. Time histories of Gp of T n at T = 2 
with improved and standard action 



following we apply Iwasaki's Improved action for 
the simulation of SU(3) gauge theory. 

From roughly 0.5 x 10 6 data, we obtain Gp for 
Tn and Tn ■ But they have still rather large er- 
rors. The fit of Gp with parameters in the spec- 
tral function given by Eq.(2) is done with SALS. 
The shear and bulk viscosities are calculated by 
these parameters by Eqs.(3) and the error are es- 



is consistent with the arguments of S. Gavin]?)] 
that in the pure gauge theory, the bulk viscos- 
ity and heat conductivity should be zero. In our 
calculation, Gp of T\\ from which the heat con- 
ductivity is calculated, has large background and 
it has been impossible to get signal from it. 

In order to know the shear viscosity in the phys- 
ical unit, we should determine the lattice spacing 
a in these /3 values. For this purpose we have 
started to determine the finite temperature tran- 
sition points p c at Nt = 8. We made a his- 
togram analysis on 16 3 x 8 lattice. The results 
for the Polyakov susceptibility are shown in Fig. 3. 

The spacial volume is too small to make a pre- 
cise determination of (3 C on our lattice, but we find 
that the transition region is 2.71 < f3 c < 2.73. 
By using the finite size scaling formula and the 
transition temperature T c ~ 276(3) (2)MeV de- 
termined by Tsukuba group Q and assuming 
asymptotic scaling for (3 > 2.73 region, the lat- 
tice spacing is estimated. 

The preliminary results for the shear viscos- 
ity in the physical unit are shown in Fig. 4. It is 
found that it increases with temperature. The T 
dependence is consistent with T 3 , because rj x a 3 
is almost independent of (3 as shown in Fig. 2 and 
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. . creases with temperarure roughly as T 3 

The bulk viscosity is consistent with zero. And 
the heat conductivity is very difficult to calculate 
because of large background for Gp of T14. 

Our results depend strongly on the ansatz of 
spectral function in Eq.(2). In order to study the 
functional form for the spectral function, we have 
started the simulation with anisotropic lattices. 

We could not calculate the Gp in the confined 
phase We think that this is because the energy 
momentum tensor operators should be written by 

1 1 the hadron fields, rather than the gluon fields. In 

order to certify this hypothsis, we have started 
the calculation of gluon propagator at finite tern- 
Figure 3. The polyakov susceptibility on 16 3 x 8. perature with improved action. The preliminary 

results on the small lattice 8 3 x 4 shows that in the 
confined phase, gauge copies are found in lattice 
version of Lorentz gauge and the gluon propaga- 
tors show the peculiar behavior, which Nakamura 
group H has found on the large lattice with stan- 
dard action. 
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Figure 4. Shear viscosity in the physical unit 



a oc 1/T when Nt is fixed. But clearly more ac- 
curate data are necessary. It is also seen that it 
seems finite around T c , which is slightly smaller 
than the typical hadron massess. What is the 
physical effects on the phenomenology of quark 
gluon plasma, when it has shear viscosity with 
this magnitudes, is a very interesting problem. 

3. Conclusions and Further Problem 

We have obtaind the transport coefficients of 
quark gluon plasma from the lattice QCD cal- 
culation by using Iwasaki's improved action, as 
shown in Figs. 2 and 4. The shear viscosity in- 
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